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ABSTRACT 

We have fit the Fundamental Plane of Ellipticals (FP) to over 400 early-type galaxies in 20 nearby clus- 
ters (cz ~ 4000 — 11000 km s _1 ), using our own photometry and spectroscopy as well as measurements 
culled from the literature. We find that the quality-of-fit, r.m.s. [log(cr)], to the average Fundamental Plane 
(FP) varies substantially among these clusters. A statistically significant gap in r.m.s. [log (a)] roughly 
separates the clusters which fit FP well from those that do not. 

In addition, these two groups of clusters show distinctly different behavior in their peculiar velocity 
(PV) distributions. Assuming galaxies are drawn from a single underlying population, cluster PV should 
not be correlated with r.m.s. [log(er)]. Instead, the clusters with below average scatter display no motion 
with respect to the cosmic microwave background (CMB) within our measurement errors (~ 250 km s _1 ), 
while clusters in the poor-fit group typically show large PVs. Furthermore, we find that all X-ray bright 
clusters in our sample fit the FP well, suggesting that early-type galaxies in the most massive, virialized 
clusters form a more uniform population than do cluster ellipticals as a whole, and that these clusters 
participate in a quiet Hubble flow. 

Subject headings: galaxies: fundamental parameters — galaxies: clusters — cosmology: cosmic microwave 
background — cosmology: distance scale 
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1.1. The Fundamental Plane of Ellipticals 

The high correlation of the structural and kinematic properties of ellipticals suggests that these galaxies formed by 
similar processes and are largely virialized. Assuming structural symmetry, isotropic velocities, (v 2 ) = a 2 , a constant 
mass-to-light ratio (M/L), and structural homology, a tight correlation is expected as a natural result of virialization : 

a 2 oc (SB)r, (1) 

where a 2 = the velocity dispersion within the galaxies, r = a fiducial radius, and (SB) the average surface brightness 
within r. In log space this expression relating basic physical properties collapses to a plane and, hence, is called the 
Fundamental Plane of Ellipticals (Djorgovski & Davis 1987; Lynden-Bell et al. 1988; Lucey, Bower, & Ellis 1991). 
In terms of observables the FP becomes 

log(r e ) = a log(er) + j3 (fi) e + Const, (2) 

with a = 2, j3 = 0.4 representing the virial plane. Standard units of measure are km s _1 for the central velocity 
dispersion, a, arcseconds for the r 1 / 4 — law half-light radius, r e (effective radius), and magnitudes per arcsec 2 for the 
average surface brightness within that radius, (/i) e (effective surface brightness). 

Empirically, ellipticals occupy a narrow range within this three-dimensional parameter space indicating they are 
a homologous family. SO galaxies also lie on the FP due to the fact that, while errors in the estimated r e and (/z) e 
are made by fitting r 1 / 4 — law profiles to galaxies with disk components, the errors cancel such that the FP quantity, 
"mboxlog(r e ) — (3 (fi) e ", is left remarkably unaffected (Saglia et al. 1997; Scodeggio et al. 1998; Kelson et al. 
1999). The observed FP, however, is misaligned with the virial prediction. Ellipticals are better described by, a ~ 1.4 
and (3 ~ 0.3. This tilt of the FP implies the assumptions used in Equation [l] are not quite true. Nonetheless, the 
uniformity of ellipticals makes these galaxies important standard candles since the observed r e provides a direct 
indication of distance (Dressier 1987; Dressier et al. 1987; Lucey & Carter 1988). The small observed scatter about 
the FP translates to ~ 19% error in the distance to individual galaxies (J0rgensen, Franx, & Kjaergaard 1996; Hudson 
et al. 1997; present work). The FP is thus comparable as a distance indicator to the Tully-Fisher (TF) relation applied 
to spirals (de Carvalho & Djorgovski 1989). 



1.2. Motivation & Goal for this Study 

The FP, therefore, should be a powerful tool for measuring relative cluster distances as well as deviations from 
the Hubble flow. However, surveys of large scale motions using the FP and TF have found significantly different 
cosmic flows. For example, the recent FP "Streaming Motions of Abell Clusters" survey (SMAC) (Hudson et al. 
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1999) reveals a large bulk peculiar motion on the galaxy cluster scale. Yet, while SMAC agrees with the 15 cluster TF 
survey (LP10K) of Willick et al. 1999, most TF surveys have shown little peculiar motion on these scales. The large 
and somewhat deeper TF cluster survey of Dale et al. 1999a limits the motion to less than 200 km s , agreeing with 
TF surveys going back to the Aaronson et al. 1986 TF survey. To further confuse the issue, Lauer & Postman 1994, 
using a novel method employing brightest cluster galaxies in 119 clusters, detected a large bulk flow of the 'Abell 
Cluster Inertial Frame" (ACIF). But the direction is in conflict with the results of the aforementioned surveys, none of 
which have the large sky coverage of the ACIF survey. 

As can be seen from the compilation of previous measurements (Table |l|), the results of these surveys are clearly 
discrepant or their errors have been underestimated. Indeed, some of the possible sources of systematic errors, as 
summarized by Jacoby et al. (1992), are still not well understood. Moreover, as emphasized by Strauss et al. 1995, 
if flows on this large a scale exist, then they can not be accommodated by any present cosmological model. These 
models assume the average random motion of clusters atop the Hubble flow is driven by gravity and the underlying 
mass distribution and thereby provides a constraint on the density parameter, ft (Borgani et al. 1997; Watkins 1997; 
Bahcall & Fan 1998). 

To address this problem, we examine the FP and rigorously test the assumption that the FP is invariant from 
cluster to cluster. If this assumption is violated, then previous surveys, which have assumed a uniform fundamental 
plane, could contain significant systematic errors. We have compiled a sample of over 400 early-type galaxy in 20 
clusters. Measured distances to clusters should be more accurate than to individual galaxies by the root number of 
galaxies observed per cluster, and thus typical distance errors to individual clusters should be of order ~ 5%. Hence, 
flows on cluster scales can in principle be accurately determined. Furthermore, this sample allows us to investigate 
the effects of possible systematic errors on a cluster by cluster basis. The rest of this paper is organized as follows: in 
§2 we briefly describe our data as a thorough presentation will be the subject of a subsequent paper; we discuss our 
FP fitting procedure in §3; results and statistical analysis are presented in §4; we finish with further discussion and 
implications of our result in §5; and summarize our conclusions in §6. 

2. The Data 

Our observations have yielded useful spectroscopy and photometry of 132 E+S0 galaxies in 8 Abell clusters. 
Spectra have been obtained with the Nessie multi fiber instrument on the KPNO Mayall 4m, and imaging with a large 
format (2048 2 ) Tek CCD on the KPNO 0.9m. Redshifts and central velocity dispersions are measured by matching the 
galaxy spectra with G and K class stellar templates, using the Fourier-quotient technique (Sargent et al. 1977; Tonry 
& Davis 1979) adapted for use in IRAF by Kriss (1992). Aperture corrections are applied to the velocity dispersions, 
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due to the fact that the fibers, which are of fixed diameter, sample larger portions of galaxies at larger distances, 
systematically lowering measured velocity dispersions. For model ellipticals, J0rgensen, Franx, & Kjaergaard 1995 
find the velocity dispersion measured within an aperture varies with aperture size as a power law of index of —0.04. 
The applied aperture correction to measured velocity dispersions is then calculated using cluster redshift as a first order 
estimate of distance and the angular diameter— distance relation with qo = 0.5. 

Photometry has been corrected for atmospheric extinction, Galactic absorption (Burstein & Heiles 1984), (1 + z) A 
cosmological dimming, and ^-correction. The ^-correction to flux scales as (1 + 1.089z), calculated for an average 
elliptical galaxy spectrum (Coleman, Wu, & Weedman 1980) through our CCD and the R filter. We work in R band 
to help minimize the effects of cluster differences in age and metallicity which can cause shifts in the FP at shorter 
wavelengths (Guzman & Lucey 1993; Gregg 1995). Effective radii are calculated by fitting an r 1 / 4 — law to isophotes, 
which have been obtained using the ELLIPSE task in IRAF. A grid of seeing corrections to r e and (/i) e is calculated 
from models convolved with the PSFs of the images. 

Merging our data with recently published observations, increases the number of galaxies to 428 in 20 clusters. 
Table || lists the literature sets as well as abbreviations used hereafter. Velocity dispersions, effective radii, and surface 
brightnesses have been normalized to the SLHS97 system by minimizing residuals of each of these parameters for 56 
galaxies in four clusters common to both samples. The details of our observations and this conversion including seeing 
and aperture corrections will be described in detail in Gibbons, Fruchter, & Bothun (2000c). For some clusters, e.g. 
A400, our sample size is limited because we exclude objects with strong emission lines and evidence of a strong disk 
component. We retain only galaxies which are clearly cluster members and which have excellent spectroscopic S/N. 
Average galaxy FP errors for clusters are comparable between the three data sets except for A400, whose spectra are 
of lower quality than the other clusters. However, excluding this cluster from our analysis does not alter our results. 



3. Fitting to All Data Simultaneously 

Once the corrections described in §2 are applied to the measured cr's and (/i) e 's, shifts in the FP should be due 
purely to the change in apparent galaxy size with distance. The FP is then fit to all the galaxy data simultaneously by 
letting cluster distance be a free parameter. Specifically, we solve a set of equations with 428 galaxies in 20 clusters, 

log(r e ) = a log(a) + (3 (fi) e + 7i , i = 1, 20 (3) 

where shifts in the intercepts, Aji, reflect the offsets in log(r e ), which translate into the relative distances between 
clusters. 

Because one is seeking high precision in relative distances it may seem most prudent to minimize the residuals 
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about the FP in the direction of the distance dependent parameter, log(r e ). This particular projection also yields the 
smallest observed scatter due to the strong correlation between r e and (fi) e . However, Lucey, Bower, & Ellis (1991) as 
well as J0rgensen, Franx, & Kjaergaard (1996) have recognized that minimizing the residuals of log(r e ) will introduce 
a bias into the fit because the errors in log(r e ) and (fj) e are correlated, since (fi) e is a function of r e (see also Akritas 
& Bershady 1996). 

Therefore, like SLHS97, we isolate the velocity dispersion, a, as the only parameter with independent errors. 
Explicitly, we write the FP solely as a function of r e , 

log((7) = f(r e ) = (log(r e ) - (n) e - 7i )/a , i = 1, 20 (4) 

and solve for the coefficients, which minimize the absolute residuals, 

£ Wr^wi (5) 

i=l,20 j=\,ni 

where rii = the number of galaxies per cluster. When the absolute residuals are minimized, galaxies which lie far 
off the main relation have effectively less weight than they would under least squares. Our use of a as the indepen- 
dent variable avoids the strong bias caused by the correlated errors in r e and (fi) e when r e is instead treated as the 
independent variable. 

We solve this system of equations allowing 22 free parameters: 20 cluster intercepts, 7/s, plus common a and j3. 
In this way, we find the average fundamental plane, FP, which best fits the entire sample of galaxies, further reducing 
the influence of anomalous galaxies. To ensure the best solution within a parameter space of such complexity (22 
dimensions) requires an efficient searching algorithm. We adopt simulated annealing, because it efficiently converges 
towards the global minimum while avoiding becoming trapped in local minima (Metropolis et al. 1953; Kirkpatrick, 
Gelatt, & Vecchi 1983; Vanderbilt & Louie 1984). The convergence is quite rapid and most of the computing time is 
spent sampling the region near best fit. 

4. Results 

4.1. The FP Derived from the Present Sample 

The best-fit FP, a — 1.37±0.04, — 0.331±0.004, is shown in Figure [l|a and is in agreement with that found 
by Hudson et al. (1997). We estimate the errors on a and (3 by fixing 7.;, varying a and f3, and constructing maps of 
chi-square versus a and j3. Our chi-square assumes the average total FP scatter (intrinsic plus measurement error) is 
well represented by the measured average r.m.s. scatter in log(er) about FP, (£}. Then 

n 

X 2 = £[log(^-FP] 2 /<£ 2 ). (6) 
3=1 
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The one sigma errors, corresponding to A\ 2 = 1, are read directly from these curves. For our fit, the average r.m.s. 
scatter in log[cj], (£) = 0.065, is equivalent to a 20% error in distance to individual galaxies. Assuming the average 
scatter applies to all clusters, errors in cluster distances range from about 2 — 8%, depending on the number of galaxies 
fit per cluster. 

However, 9 of the 20 clusters, about half the total sample, show significantly more scatter than the other 1 1 
(hereafter £ denotes FP r.m.s. scatter). When these clusters are fit separately, they also fail to define a tight FP. We 
must therefore investigate the significance of the larger scatter: i.e., whether or not the high-£ fits are truly poor. The 
degree to which the FP describes the individual clusters can be seen in Figure ^p. Assuming galaxies in a cluster 
are chosen randomly from a single population, we calculate for each cluster the significance, or probability, P, of a 
chi-square this large about the FP and plot this value against £ Immediately evident is a break in P coincident with a 
separation in £ about (£). While the clusters should uniformly fill probability space, a gap this large anywhere in P 
has less than a 0.1% chance of occurring (Similarly, along the £ axis, one would expect the points to be clustered about 
(£) rather than showing the deficit that is observed). However, as this estimate of probability relies on the assumption 
that the FP is gaussian, we have also estimated the probability of this gap by simulating many trial sets of clusters from 
the total sample of observed galaxies and calculating the distribution of the maximum gap. Doing so again implies the 
observed situation is improbable, with an expected occurrence of a gap this size of 0.2%. This strongly suggests that 
our sample has not been drawn from a single galaxy distribution. 

We therefore refit for FP based on the 1 1 best-fit clusters (Figure ^a). Although the FP coefficients are only 
slightly different, a = 1.39±0.04, = 0.335±0.005, and (£) changes by - 10%, the 9 clusters with £ > (£) 
clearly appear to be outliers, falling in the upper 5% tail in probability space (Figure ||b). If truly deviant, these 
clusters may not be reliable objects for FP distance work. But do the high-£ clusters distinguish themselves from the 
rest in any other way? PV is a second, independent dimension in which we can examine the behavior of this sample 
and is after all the quantity we aim to measure. 

4.2. Measuring Distances and Peculiar Velocities 

To derive peculiar velocities, we must first fix the scale between FP relative distances and absolute distances. 
That is, we must set the conversion between 7 and cluster redshift. We begin with the assumption that the clusters are 
at rest with respect to the Hubble flow, i.e. at rest in the CMB frame. Cluster line-of-sight velocities are translated 
from heliocentric redshifts, as determined from our spectra, to redshifts in the CMB frame using a vector in the 
direction I = 264?4, b = 48?4 in Galactic coordinates with a magnitude of 369.5 km s~\ (Kogut et al. 1993). Taking 
into account (the small) cosmological corrections, we then find the angular diameter — distance zero-point which 
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minimizes the one-dimensional r.m.s. peculiar velocity, <7i£>, of the clusters. This normalization does not change 
significantly if we minimize a^u for the eleven best-fit clusters instead of for all twenty. Interestingly, this procedure 
places the Coma cluster nearly at rest, so we would have recovered a similar result had we anchored the physical scale 
by assuming no peculiar motion for Coma (which has traditionally been done). 

Finally, plotting PV against E, we uncover an astonishing result (Figure |J). There is an obvious discontinuity in 
the scatter in cluster PV which occurs across the same region of E-space as the gap in probability already discussed. 
The clusters which fit FP well are at rest with respect to the CMB with reduced x 2 < 1 while the poor-fit clusters with 
£ > (£} show large scatter in their PVs, \ 2 = 4.4. |^| Results from both the 20 cluster FP and 1 1 cluster FP are similar 
as can be seen in Figure f|. 

While one might naively expect a correlation between S and the measurement error of PV, this is not the case. If 
cluster ellipticals do indeed form one population, then the sample mean and sample standard deviation of a randomly 
drawn subset will be uncorrelated. Therefore, PV should not reflect the value of E. Similarly PV errors do not 
correlate with S; we therefore scale the errors by (£}. As we will show in the following section, this null hypothesis 
is heavily disfavored as the observed non-uniformity in g^u is highly significant. 



4.3. Significant Variation in 0\d 

To examine the significance of the variation in a\o, we have used two approaches, linear regression and analysis 
of variance (ANOVA). Simple regression favors a strong correlation between PV and Ewhether we use (£) to estimate 
the PV errors (null hypothesis), or use the individual S. Straightforward weighted regressions of "|PV| on S" and "E 
on |PV|" yield slopes which are nonzero with 99.9% confidence. The true slope lies somewhere between the forward 
and inverse regressions. But because of the two-way errors, it is difficult to decide significance analytically, so the 
problem can also be dealt with numerically. Resampling by complete permutation along the E axis is a measure of the 
likelihood of the observed correlation. Permutation ignores errors in E and complete mixing of the 20 points without 
replacement is allowed. Under the assumption of no correlation, the expected slope is zero, negative and positive 
slopes being equally likely. The resultant distribution shows the observed slope to be unlikely, with a probability of 
less than one percent. [] 

The null hypothesis likewise fails ANOVA tests. The measured o\d for the two subsamples of clusters (defined 

3 Fractional errors in distance are predicted based on {£} and the FP coefficient, a, and scale inversely with y/n ga i. Estimated PV errors are 
distance errors added in quadrature with redshift errors. 

4 Note, however, if two separate populations do not show such a correlation, this test is inappropriate, in which case ANOVA is more meaningful. 
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by the gap in S) are not consistent with a single PV distribution. The ratio of chi-squares disfavors the two groups 
being drawn from the same distribution with probability 99.94%. [] Again permuting the points along E— space allows 
us to estimate the probability of the present result without an assumption about the nature of the distribution. Using 
this bootstrapping procedure one finds the significance for the ratio of observed chi-squares is again above the 99 th 
percentile. The behavior of the sample both in S— space as well as in PV— space tends to contradict at very high levels 
of significance the assumption that we are observing a single population of clusters. Taken together, these statistics 
are convincing evidence that we are seeing two populations. If this bimodality is an intrinsic property of clusters of 
galaxies, then all current FP samples potentially mix together these two apparently different populations. 

5. Errors in FP Distances or Peculiar Velocities? 

We have established that our sample of clusters divides into two groups defined by significantly different intrinsic 
scatter about the FP. Additionally, we find that these two groups have radically different PV distributions. The obser- 
vation that only the high-S clusters show peculiar motions implies that such clusters can inject significant error into 
the measure of a\r>. This could have the consequence of producing a spuriously large PV signal for the entire sample. 
Alternatively, these two groups could indeed be subject to truly different kinematics. Either scenario, suggests that the 
nature of the inferred PV field is sample dependent. 

We find, however, no evidence that this apparent separation of the clusters into two populations is due to system- 
atic problems with our data. There is no correlation of S with cluster distance (e.g. more distant clusters do not have 
poorer fits). We see no effect due to data quality or data source; the clusters which we observed ourselves and the 
sample which was culled from the literature are equally divided between the two £ groups, and with the exception of 
A400, the quality of the imaging and spectra of all the clusters are comparable. We are also confident that our cluster 
samples are not significantly contaminated by field ellipticals or SOs; we are sampling the centers of clusters over an 
area of approximately one half square degree and our redshift criterion is \z—z mean | < 3cr, where a — cluster velocity 
dispersion. We calculate that, with the possible exception of 7S21, the expected number of field interlopers per cluster 
in our sample is far less than one. Addition of an Mg2 line index term to diminish cluster to cluster differences in the 
age/metallicity of the stellar populations also does not improve the fits. We find no strong evidence for curvature along 
the FP. In addition, the galaxies in each cluster cover similar ranges in log [a) so that we expect any error introduced 
by sampling different portions of the FP is at most a minor effect in our final fits and distance measurements. The 

5 This is not the traditional F-test because the comparison is not between two independent fits. For this problem, the number of degrees of 
freedom for the subgroups is us"t 1 x nsu bgroup. Specifically for nclust = 20, the d.o.f. are 8.55 and 10.45 for nsubgroup = 9 and 11 
respectively. Furthermore, the ratio of \ 2 is no longer purely F-distributed, so we have computed the proper cumulative distribution function. 
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high-E clusters do on average have a smaller number of observed members than low-E clusters. However, we would 
expect a small sample size to only increase the scatter of E and not bias its value. Bootstrap tests on subsamples of 
galaxies chosen at random from Coma show this to be the case. Rather, the high-E clusters are also on average the 
least rich clusters, and as we discuss below, there is good reason to believe that E is a function of cluster richness. 

If our study is relatively free from systematic errors, then our statistical analysis shows that different clusters 
have different amounts of intrinsic log(cr) scatter about the FP. But what property of a cluster might be driving this 
observed difference? We have looked at many cluster properties as predictors of whether a cluster will be a high- or 
low-E cluster and have found only one strong indicator : cluster X-ray luminosity, as can clearly be seen in Figure || 
The four brightest X-ray clusters, those with L x > 10 44 erg s _1 , all have low E. This result suggests that early-type 
galaxies within the most massive, well-virialized clusters form a more homogeneous population than cluster ellipticals 
as a whole. 

A weaker, but similar correlation, can be seen with cluster spiral fraction, which is also an indirect indicator of 
the dynamical evolutionary state of a cluster (see for example Dressier et al. 1997). Clusters with E > 0.08 exhibit a 
spiral fraction of 55 ± 1% compared to 42 ± 3%. This trend of E with spiral fraction is shown in Figure ^. Similarly, E 
appears to be related to the cluster velocity dispersion, which is itself, of course, closely tied to cluster X-ray luminosity 
(Figure With one exception, clusters with velocity dispersions less than 400 km s _1 tend to have high E whereas 
those with dispersions greater than 800 km s _1 tend to have low E. 

There are also indications that substructure may be an important component of the increased scatter observed 
in some clusters. Take the case of the relatively spiral rich, low velocity dispersion clusters A2151 and A400. In 
A2151, the early type galaxies in the core have a mean velocity offset of —900 km s _1 from that of the spirals 
(Zabludoff et al. 1993, Dale et al. 1999b; this paper). Indeed, A2151 and A2147 are clearly entangled in a mutual 
interaction (e.g. Bird et al. 1995; Maccagni et al. 1995), with a component of peculiar velocity reflecting infall and/or 
merging of subclusters. A similar offset between mean elliptical and cluster velocity has been observed in A400 
by Beers et al. (1992). In addition, both A400 and A2151 contain many examples of interacting/peculiar galaxies 
(see Bothun & Schommer 1982), indicating the prevalence of low velocity dispersion subgroups within the overall 
structure. Furthermore, substructures that are separated by ~ 10 Mpc or less will not be adequately resolved in 
kinematic distance space by the current relative indicators. An increase in E can be introduced due to a range in 
distance along the line of site. Indeed, such a spread in distances has already been proposed for one rich, but high E 
cluster, A2634 (see West & Bothun 1990; Scodeggio et al. 1995). 

One reason for the fact that highly X-ray luminous clusters do not show significant peculiar velocity may be 
that X-ray luminosity coincides with the virialization of infalling structures. Whatever the source of the consistent 
behavior of these clusters, it appears that X-ray luminosity provides a good pre-filter for the selection of clusters to be 
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used in PV studies. As will be shown in paper II of this series (Gibbons, Fruchter, & Bothun 2000b), such a pre-filter 
produces a sample that exhibits no large scale flow. 

6. Conclusions 

Using high-quality data and a robust fitting technique, we have determined that clusters of galaxies show a large 
range in quality-of-fit of their member galaxies to the FP. Statistical tests described in §4 argue that the observed 
separation in £ is significant at the ~ 10 -3 level. Additionally, an examination of the distribution of clusters along the 
PV axis provides further direct support for this separation. If clusters are truly random realizations of galaxies drawn 
from a single population, then the total sample should have, on average, a uniform am- Instead, when the sample is 
divided on the basis of the observed scatter about the FP, the two subsamples have drastically different properties in 
the independent dimension provided by the PV. Again, this effect is statistically strong at more than the 99% level. 
Simply put, high-E clusters have large PV and low-E clusters have essentially no PV, within the observational errors. 
Taken together, these two results cannot be reconciled with the idea that these clusters are randomly drawn from a 
single population. 

The behavior of E appears to be an intrinsic cluster property. We have found that a subset of the clusters, the 
X-ray luminous clusters, exhibit a well-defined FP. This suggests that the most massive and virialized clusters possess 
a homogeneous population of ellipticals and therefore provide the most accurate measures of relative distance. About 
half of the cluster sample possess £ below the mean. When using this half of the sample, we find the clusters to be 
at rest within the CMB frame (see details in the forthcoming Gibbons et al. 2000b). The high-E clusters provide the 
bulk of the positive PV signal when all the clusters are averaged together. 

Overall, our results strongly suggest that making a cut on the r.m.s. scatter about the FP, or more strictly on 
cluster X-ray luminosity, will lead to a more reliable sample from which bulk flow properties can be established. The 
current level of disagreement between various observations of large scale flows (e.g. Table [j]) may possibly reflect the 
fact that less reliable samples have been used. 
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Survey 


Method 


PV w.r.t. CMB (km s" 




Survey Depth (km s 


a id (km s 


ACIF 


BCG 


689 ± 178 {I = 343 b = 


52) 


15,000 


- 400 


SMAC 


FP 


630 ± 200 (I = 260 b = 


-1) 


12,000 


- 350 


LP10K 


IRTF 


720 ± 280 (I = 266 b = 


19) 


12,000 


- 415 


Dale et al. 1999 a&b 


IRTF 


< 200 




18,000 


341 ± 93 


Watkins 1997 


IRTF 






12,000 


265+^ 6 



Table 1: Recent Surveys : These surveys disagree in both the magnitude and the direction of their derived flows. 
References appear in the text. N.b. The Watkins analysis is of subsets of the Giovanelli et al. (1997) and Willick et al. 
(1995) surveys. 
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1 


0.088 
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J8 


150.3 
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1 = this work; 2 = Lucey & Carter 1988 (LC88); 

3 = J0rgensen, Franx & Kjsergaard 1995a,1995b (JFK95a,JFK95b); 

4 = Lucey, Guzman, Smith & Carter 1997 (LGSC97); 

5 = Smith, Lucey, Hudson & Steel 1997 (SLHS97) 



Table 2: The Cluster Sample : column (1) cluster name; (2-3) cluster position in Galactic coordinates; (4) heliocentric 
redshift; (5) cluster velocity dispersion; (6) number of observed galaxies; (7) data references as defined above; (7) 
scatter about the average FP; (8-9) FP parameters for the fits to individual clusters. 
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Fig. 1. — (a) The FP for the entire sample of 428 galaxies. This projection shows the scatter about log(r e ). (b) The 
significance of x 2 about FP for each cluster vs. FP scatter. The gap seen in probability space is highly significant, 
making it unlikely that all clusters have been drawn from a single population. 



-16- 




Fig. 2. — (a) The FP for the 1 1 best-fit clusters, (b) Same as Figure [l]b but for the 1 1 cluster FP fit. When fit this 
way the high-S clusters are outliers, falling in the upper 5% tail in probability space, while the rest are consistent with 
being uniformly distributed. 
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Fig. 3. — Cluster PV vs. S. Each point represents a cluster containing ^10 — 80 observed galaxies. The x-axis is the 
scatter within the cluster with respect to FP. The abscissa is the PV as derived from redshift minus FP distance. The 
reduced x 2 about zero PV for the entire sample is shown at the top of the figure. (£) (vertical line) roughly marks 
the separation of the sample into two subgroups, over which \ 2 rises dramatically (values shown at the bottom of the 
panel). The errors on £ are shown in the bottom panel for clarity. One can see a paucity of clusters about (£), where 
one would expect the density of clusters to be highest if clusters are drawn from a single population. 
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Fig. 4. — Comparison of the PV distribution for the 20 and 1 1 cluster fits. As is readily seen, use of one fit or the other 
does not significantly change our results. 
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Fig. 5. — S vs. X-ray luminosity, L x . Note that the massive, virialized clusters fit the FP well, while the quality-of-fit 
for the lower mass clusters cannot be predicted. The most X-ray bright clusters (Lx <; 10 44 erg s^ 1 ) are labeled. 
The clusters for which there exist no X-ray data are shown along the left side of the panel. Lx from White, Jones, & 
Forman (1997). 
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Fig. 6. — S vs. cluster spiral fraction (SF). A weak correlation between these two variables is evident and is in the 
same sense as expected from the X-ray data. Spiral fractions from Dressier 1980. 
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Fig. 7. — E vs. cluster velocity dispersion. The high dispersion (;> 800 km s 1 ) clusters have low E. 



